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Abstract 

We predict the dependence on energy of photo (electro) production processes: 7(7*) +p — > V + X 
with large rapidity gap at small x and large momentum —t transferred to V in pQCD. Here V is a 
heavy quarkonium (J/ip, T) or longitudinally polarized light vector meson (in the electroproduction 
processes), etc. In the kinematics of HERA we calculate the dependence on energy of cross sections 
of these processes as the function of momentum transfer t, photon virtuality Q 2 and/or quarkonium 
mass. In the kinematical region Qq < — t <C Q 2 + My the nontrivial energy dependence of the cross 
section for the vector meson production due to the photon scattering off a parton follows within 
QCD from the summing of the double logarithmic terms. In the second regime -t > Q 2 +M v within 
DGLAP approximation in all orders of perturbation theory the qq — parton elastic cross section 
is energy independent. We show that the correct account of the double logarithmic terms and of 
the gluon radiation including kinematical constraints removes the disagreement between pQCD 
calculations and recent HERA experimental data. The explicit formula for the dependence of the 

differential cross section—— — of these processes on s 7 *at is obtained. We show that perturbative 

LLLLLJu J 

Pomeron type behavior may reveal itself only at energies significantly larger than those available 
at HERA. In addition we evaluate the energy dependence of DVCS processes. 
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I. INTRODUCTION. 



The hard inelastic exclusive photo(electro) production processes: 7*p — > V + 
rapidity gap + X with a large rapidity gap and proton dissociation into hadrons draw a 
lot of attention recently (V is vector p— meson, charmonium J/ip, or bottomium) . (Such 
processes can be dubbed hard inelastic diffractive (HID) processes, in order to distinguish 
them from the diffractive processes, when the final particle remain a proton: X = p). These 
processes with V = J/ip were extensively studied experimentally at HERA [jJli] due to a 
reasonably large cross section and a clean final state - production of an isolated charmonium 
state with large transverse momentum which is separated from the proton remnants by a 
large rapidity gap. Theoretically this process is advantageous, since the large charmonium 
mass my makes it possible to evaluate cross section of the process in perturbative QCD, 
for different values of transverse momentum transfer —t and virtuality of the photon Q 2 
0-01 • Moreover, the QCD factorization theorem has been proved 0, ls[ for diffractive photo 
(electro) production of vector mesons at least for the case of the longitudinal photon putting 
the calculations of diffractive cross sections on the solid theoretical basis. Processes where 
proton remains intact are strongly suppressed at large t by the nucleon form factor and will 
be considered briefly only. 

Moreover, the charmonium photoproduction is to be probed soon at ultraperipheral col- 
lisions at LHC |9j for significantly larger energies. The CMS and ATLAS detectors are 
well suited for such an investigation since they cover large rapidity intervals. The ALICE 
detector may be capable of studying this process in a certain rapidity range as well. 

Phenomenological analysis of the HI and ZEUS data on diffractive charmonium pro- 
duction [l], 0], carried out in ref. [Io[ within framework of the QCD factorization the- 
orem has demonstrated that the experimental data can be described if the energy de- 
pendence of two body pQCD scattering amplitude would be significantly slower at large 
—t, than at t = 0. Recently, a new experimental data on the photoproduction process 
7P — > J/ip + rapidity gap + X was reported [3j . This data also indicate a rather dramatic 
slowing down of the dependence of the cross section of this process on the energy with in- 
crease of —t. This conclusion is consistent within the experimental uncertainties with the 
phenomenological analysis of ref. [Io[ • The observed behavior differs 0] from the predictions 
which assumed that the energy dependence of the amplitudes describing two body processes 
f(s, M\, t, Q 2 ) can be evaluated in terms of the BFKL approximation with the coupling 
constant independent oft 0], or with a e // = a s (i) (llj . 

The novel QCD effect discovered in this paper is that the pattern of energy dependence of 
cross section of the HID on t, observed at HERA, follows from the DGLAP approximation 
of pQCD. The amplitude of the hard two body collisions has a rather simple form in a wide 
interval of high energies where the double logarithmic (DL) approximation is legitimate: 

/ = fDL(x u x 2 , M 2 x ,t, Q 2 , M 2 )F(x, M x , t). (1.1) 



Within the DGLAP approximation [12Hl4j the function foi is obtained by summing large 
d s {N c /2'w) ln(xj/x) ln(Q 2 + M v )/(Ql — t) terms, that arise in the integration over parton 
transverse momenta in the domain where a s < 1. Here xj = —t/M x is the fraction of the 
proton momentum carried by the parton involved in the large t elastic scattering. Variables 
Xi are the fractions of the proton momentum carried by gluons exchanged in t channel which 
are attached to a parton of the proton (For the explicit definition of xi,x 2 see section 2.). 
Since at large t transverse momenta of exchanged gluons are large, Xi are not vastly different. 
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Note that it has been demonstrated in ref. 15], see also refs. 13, T^, that increasing with 
energy DL terms provide a good description of the structure functions of proton measured at 
small x at HERA. This observation suggests that at HERA energies F = 1. ( We normalise 
foL an d F in such a way that F = constant = 1 at low energies.) 

The function F = 1 within the LO DGLAP approximation to pQCD, but it is a function of 



energy within the BFKL and resummation models 17-19]. This is because these approaches 
take into account the contributions that increase with the energy but do not contain double 
logarithmic (DL) terms. We estimate energy dependence of F in the paper. The effect is 
small at HERA energies but may become noticeable at significantly higher energies. 

In the kinematic region Ql < — t < Q 2 + My the energy dependence of hard amplitudes 
is determined by foL- The equation 11.11 differs in this kinematic range from the Pomeron 
exchange expression (s / ' so) ap< ^ } ~ 1 often used to describe the data. In particular, DL ap- 
proximation predicts strong dependence of phenomeno logical ap(t) on Q 2 and on incident 
energy (i.e. non universality, dependence on the external conditions). This non universality 
of ajp specific for DGLAP approximation has been observed at HERA at t — |20j: 

F 2 (x,Q 2 ) ~exp(r(Q 2 )log(l/a;)),r(Q 2 ) ~ 0.05 log(Q 2 /A 2 ). (1.2) 

In the kinematical region — t > Q 2 + M v , Jjjl is energy independent. Hence the energy 
dependence of the cross section is determined by the function F which does not contain 
large logarithms from the integration over parton transverse momenta and which is equal 
to one at HERA energies. Thus pQCD predicts a sharp decrease with —t of the rate of 
the rise with energy of the photo(electro) production cross sections as compared to forward 
scattering at HERA energies. 

The estimates of the kinematic range dominated by double logarithms [l8|, 2l|| show 



that the universal Pomeron behavior may be valid only for the energies well beyond the 
kinematical region: ln(xj/x) ~ (2 3) \og(Q 2 /(— t + Qq)) occupied by double logarithms, 
i.e. when the DL terms disappear. At Q 2 + M v > lOGeV 2 this condition corresponds to the 
kinematics far beyond the kinematical range of HERA. The explicit analysis of the phase- 
space constraints on the multi Regge kinematics due to the energy-momentum conservation 
shows that even at ultraperipheral processes at the LHC (where x down to 10~ 6 can be 
reached) these constraints limit the possibility of the onset of the Pomeron behavior to the 
kinematic range — t > Q 2 + M v . It will be very interesting to look for onset of a Pomeron 
behavior in pQCD in ultraperipheral collisions in this particular kinematic region. 

Let us note that the dominance of the DL terms at small x in a wide kinematical range 
is a common feature of DGLAP, one of saddle points in the improved BFKL, and the 
resummation models (see the discussion in refs. [HI, [H, [HI). In the improved BFKL 



approach [23| double logarithmic terms are accounted properly in the same way as in the 
resummation models. Thus the energy dependence given byfoi should disappear with 
increase of —t in all approaches as long as the contribution of the BFKL saddle point can 
be neglected. In other words, our results will remain quantitatively correct also in these 
approaches if one treats correctly the collinear terms and includes the energy momentum 
conservation constraints. The difference may reveal itself only at the energies far beyond 
HERA. On the contrary the contribution of the BFKL saddle point [l7[ in the inverse Mellin 
transformation has no DL terms and therefore it does not lead to a peculiar dependence on 
t of the rate of increase of the cross section with energy discussed in the paper. 

We solve the DGLAP evolution equations in the form of the DL approximation in a large 
interval of t. It is known that the DL approximation gives quantitative description of the 
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HERA data on the structure functions of a proton 15j . Technically, the main effect related 
to the increase of —t is the following. In the case of t = the strong energy dependence of 
QCD cross section arises because the energy dependent terms \og(xo/x) are enhanced by the 
factor log(log(Q 2 /A 2 )/ log(Qo/^- 2 )) where Qq is an initial virtuality. This factor originates 
from the integration over transverse momenta in the range Qq <C k 2 <C Q 2 . For —t 3> Qq 
the effective range of integration changes to — t <C k 2 <C Q 2 + M v , leading to the decrease 
of the rate of the energy increase with the increase of —t. At — t ~ Q 2 + M v the DL terms 
disappear all together. Then the energy dependence may arise due to either nonperturbative 
initial condition for QCD evolution or the terms that do not contain integration over large 
transverse momenta. 

Since the intermediate state invariant mass of the cc system at large —t is significantly 
larger then the masses of charmonium states we can approximate / at t ^ with unequal 
Xi by the nonforward parton distribution with x = (x\ + x<i)j1. It was shown that at small 
x and large scale controlling hardness of the process this approximation does not introduce 
significant uncertainties [tI liil. 25 



The nonforward parton distributions for diffractive Z boson production were studied 



also in ref. 26| in the kinematical region: log(Q 2 /Qo) > log(xo/x) which differs from 
the kinematical region of HID (\og(Q 2 /Q 2 ) ) <C log(xo/x)) considered in this paper. Cross 
section of this process is strongly suppressed at large t as compared to that for HID by 
the square of nucleon form factor. Moreover impulse approximation=QCD factorization 
should be modified at large t where two gluon scattering off two different partons becomes 
preferable. Besides, the dependence on the running coupling constant derived in our paper 
within DGLAP approximation, differs from the one suggested in ref. [26| . 

The onset of the black disc limit for the processes with large rapidity gap may occur in the 
kinematics where xj = —tjM\ would be sufficiently small, i.e. in the kinematics where the 
gap between V and the system X is relatively small. The distinctive signature of this regime 
is a significantly slower decrease of cross section with —t as compared to pQCD regime. Note 
also that in the black regime inelastic diffraction is suppressed. In the discussed process this 
will lead to a further slowing down of the increase of the cross section with energy at fixed 
xj,t or perhaps even to a decrease of the cross section with energy. This effect is amplified 
with increase of xj since the transverse distribution of gluons becomes more localized in the 
transverse plane. 

The paper is organized in the following way. In Sec. II we discuss the kinematics of the 
HID process. In Sec. Ill we consider the DGLAP evolution equations for the nonforward 
GPD and solve them in the double logarithmic approximation. In Sec. IV we evaluate 
the energy dependence of the amplitudes of hard diffractive processes for the kinematics 
achieved at HERA . In Sec. V we discuss briefly dependence on energy of the cross section 
of DVCS process as a function of t. In the conclusion we discuss the directions for the future 
progress. 



II. KINEMATICS OF THE DIFFRACTIVE PRODUCTION AT HIGH ENER- 
GIES. 

The theoretical framework for the calculation of inelastic diffractive processes in the high 
energy processes has been developed in the seventies and it is the so called triple-Pomeron 



limit [27]]. The original approach used independence of Pomeron tr aject ory on external 



conditions - the property of the soft QCD. It has been understood later [28] how these ideas 
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can be adjusted to evaluate the cross sections of the hard processes with large rapidity gap 
directly in pQCD where the QCD factorization substitutes the Pomeron pole factorization. 
The analogue of the triple "Pomeron" vertex is calculable within the pQCD (see Fig. 1 
and discussion below). The natural variables in these processes are the the invariant mass 
of hadronic states produced in the proton dissociation Mx, the square of the transverse 
momentum transfer —t = — (p 7 — py) 2 and the square of the invariant energy for qq parton j 
elastic scattering: 

S = Sqq-\-" parton j" Q > (^"-Q 

s = W^ p is an invariant energy squared of a full proton-photon system. These quantities are 
connected by the kinematic relation 

xj = -t/(M 2 x -m 2 p -t), (2.2) 

The rapidity gap for sufficiently large momentum transfer -t and invariant energy s is given 

by 

5y = In— == S . (2.3) 

y/(M$-t)(M* x -t) 

Finally, there exists a kinematical boundary on —t. For the forward scattering in the kine- 
matics where M\ 3> Q 2 : 

-t min = (M 2 x -m 2 p )(M 2 +Q 2 )/2s. (2.4) 

This boundary is however irrelevant from the practical point of view, since in the kinematics 
characteristic for the processes with large gap in rapidity M x /s <C 1. Therefore in the 
essential kinematical domain investigated in this paper: \t\ 3> \t\mm an d corrections due to 
t m i n can be neglected. 

Note that though in principle Mx can be measured using energy-momentum conservation 
in terms of the momentum carried by vector meson, in practice it can be determined from 
the measurement of the rapidity interval occupied by system X. 

Cross section of the HID processes — - — is calculable at large —t within pQCD as the 

dtdxj 

consequence of the QCD factorization theorem cf. [4j, |5j, |lOj : 



da da, 



dtdxj dt 



(81/16)G(xj, t) + £te(xj, t) + q t (xj, t)) . (2.5) 



The cross section is the product of two factors, both of which dependent on W 1P : G(xj,t) 
and the cross section of photon scattering on a parton j given by a nonforward amplitude 

/ (see Fig. 1) (Here 7+ ^ y+J is the cross section of scattering off a parton j ). In order 

to evaluate the energy dependence of the cross section one needs to take into account the 
energy dependence of both factors in eq. 12.51 
The first factor in eq. 12.51 is 

da \f\ 2 , s 

— = — . 2.6 
dt 16tt v ; 

Here the amplitude s'f is a hard amplitude of a photo/electro production of a system V 
when a virtual photon scatters off a parton with 4-momentum xjp. The amplitude f is 
a convolution of impact factor describing transition 7* — > V with the generalized parton 
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distribution D where the initial condition for the QCD evolution is the amplitude of the 
scattering off a single parton. This amplitude depends on four parameters —t, Q 2 , M v 
and the effective energy of parton-photon system (eq. 12.11) . In addition D depends on two 
arguments Xi rather than on xj\ 

Xl = (M 2 + Q 2 )/(s' + Q 2 ); x 2 = (M 2 - M v )/(s + Q 2 ). (2.7) 

Here M 2 is the invariant mass of intermediate state in the impact factor. In the charmonium 
production M 2 is approximately the invariant mass of cc pair: 

M 2 = k l±^ (2.8) 
z 1 — z) 



c 



and z is the fraction of the photon momentum carried by one of the charmed quarks, m, 
is the running charmed quark mass. k t is typical momentum of the charmed quark. To 
account for nonzero t = —A 2 one should substitute in the above formulae k t by k t — zA. As 
the consequence of large mass of c quark the contribution of the end point (z ~ 0, z ~ 1) is 
negligible. So essential z w 1/2. 

Since M 2 is significantly larger than M v for large A 2 , Xi are not very different. So it can 
be demonstrated following ref. 0, [25| that to a good approximation D is equal to diagonal 
parton distribution, but with x = (xi + x 2 )/2: 

x l D(x 1 , x 2 , t, Q 2 , Ql) « [(an + x 2 )/2\D{{x 1 + x 2 )/2, t, Q 2 , Q 2 ). (2.9) 

x ~ (2M 2 — M v + Q 2 )/2xjs. (2.10) 

Technically this result is due to the fact that the initial condition for the generalized parton 
distribution for small x practically coincides with a diagonal one with x = {x\ + x 2 )/2 [25]. 
The main contribution to skewness comes from the cell in the Feynman diagram of Fig. 
1 closest to the vector meson. It is easy to check that this is a property of the DGLAP 
dynamics in the kinematics -t < Q 2 + M v . For — t > Q 2 + My there are no logarithms 
from the integration over transverse momenta and nonforward distributions discussed below 
do not depend on energy in the DGLAP approximation. Thus the energy dependence of / 
and nonforward diagonal distribution are the same. 

The study of the first factor in eq. 12.51 will be the central subject of our paper. The total 
cross section is given by the integral 

dt Jkinematical cuts CbtcLXj 

We carry on our calculations for photoproduction in the kinematic region — t min <C — t < 
My and for electroproduction in the kinematic range — t m[n <C — t M v + Q 2 . We shall 
see below that for complimentary kinematic range of —t > My + Q 2 the energy dependence 
is quite different. 



III. DGLAP EVOLUTION EQUATIONS FOR NONZERO -t. 
A. General structure. 

Let us now consider the nonforward parton distribution that enters in eq. 12.51 In this 
section we calculate it in the LO DGLAP approximation and explain how to generalize these 
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results to the NLO approximation of pQCD. It was explained above that the generalized 
parton distribution can be approximated at small x by the parton distribution with the 
argument x = {x\ + x<z)/2 and with the transverse momentum transfer t = —A 2 . We now 
proceed to calculate this distribution. First, it is easy to show that in the class of gauges 



(CA) = 0, where C^ = Cip^ + c-i(f (q are numbers of the order 1), used in ref. [13|, [16| only 
ladder diagrams depicted in Fig. 2 contribute to the parton distribution. (This is not true 
in other gauges, in particular in the Feynman gauge). To achieve similarity with the parton 
model description of this ladder we will take in the paper C2 = 0. Throughout this section 
we will consider the partons in the ladder to be massless and Q 2 is parameter characterizing 
hardness of the process: photon virtuality and/or large mass of quarkonium. 

The calculation of the nonforward DGLAP ladder goes in the same way as for forward 
DGLAP ladder. In order to find the contribution of a given ladder cell i we use the Sudakov 
variables: 

hi = a { q + fop + k\ (3.1) 
The gluon propagator in the gauge CA = is given by 



^ ' k 2 + i0 



^ cw + c v w c 2 

9flu Ck + JkCf 



(3.2) 



This gauge is known to be free of ghosts. Moreover, it was shown in ref. [16] that in the 
leading logarithmic approximation (LLA) the contribution of the term k^k v is zero. As a 
result, since we are working in a LLA we can use the propagator 13.31 



D*(k) = 5 ab 



k 2 + i0 



_ C»k v + C V W 
9 » v Ck 



(3.3) 



The proof [16| is valid for —t min <C —t <S Q ' + My In the kinematical region s ^> — t ^> 
Q 2 + M v this prove should be valid as well since the contribution of the new structures oc 
is suppressed by the power of energy. 

We now perform the standard algebraic calculations using the propagator 13.31 for the 



cell i. In the DGLAP approximation k 2% ^> and the integrals over longitudinal 

and transverse momenta decouple. The use of the gauge invariance allows to demonstrate 
that only one tensor structure, g* = g^ v — (p u q fl + p^q v )/pq leads to the contribution 
containing terms oc \a(Q 2 /Ql). There are additional tensor structures that appear for the 
nonforward ladder but they do not lead to the logarithmic contributions in the integral over 
the transverse momenta. Moreover, integrals over transverse momenta of gluons exchanged 
between different cells do not produce \a(Q 2 /Ql). Thus the same structure leads to large 
logarithms in both forward (t = 0) and nonforward cases. 

Moreover, the direct calculation shows that the t dependence is present in the integrand 
only in transverse momenta of propagators in the cell. Hence the DGLAP kernels $(z) do 
not depend on t. The calculation of integral over «j is done by taking residues. As a result 
the contribution of a single cell i in the DGLAP ladder can be written in the form 

dP t = a s (Q 2 )dz$(z)I(A)gl u . (3.4) 

Here Z{ = fy/ In deriving eq J3.4l we took into account the gauge invariance, and summed 
the s and u channel contributions (that gives a factor of 2). The integral over transverse 
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momenta /(A) has the form: 



/(A) = Lj d% 



1 

+ 7^ 



= log(Q 2 )/(<5o ~~ ^)) + non logarithmic terms, 

(3.5) 

where t = —A 2 . The integral over k 2 is over the range Q\ to Q 2 . However the presence of 
large —t leads to a cancellation of the contribution of small k 2 < —t. Hence the contribution 
of a single cell i to the ladder has the form: 

dP t = a s dz<S>(z) \og(Q 2 / - t + Q 2 )g%, (3.6) 

instead of 

dP % = a s dz$(z) log(Q 2 /Q 2 )g t lMU , (3.7) 

for the DGLAP ladder with t = 0. Here &(z) is the conventional DGLAP kernel. The full 
answer for the ladder is the convolution: 

E( a s log(Q7(Qg-*)) rIL [<f m) .„ (3 .8) 

The main difference (besides the definition of x -cf. Eq. 12 . 10}) between nonforward and for- 
ward distributions is the change of the argument of the logarithm arising from the integration 
over transverse momenta. 

Note that the latter expressions were derived in the kinematics Q 2 ^> — t ^> Qq, while in 
the kinematical region Q 2 ~ —t one cannot single out the logarithmic contributions since 
they cannot be considered any longer as a large parameter. This means that the upper limits 
of integration in the integral of eq. 13.51 must be taken in such a way that for —t ~ Q 2 the 
argument of the logarithm is equal to one. Such choice is allowed since we work in the LLA 
and all integrals are calculated with logarithmic accuracy. In addition, the expression for the 
nonforward distribution should smoothly match with the forward distribution as — t — > 0. 
Thus the arguments of the logarithms that appear in the nonforward DGLAP distribution 
are indeed \og(Q 2 /(Ql — t)). 

We calculated only the imaginary part of the nonforward ladder. The real part is also 
nonzero and can be reconstructed using the dispersion relation over the energy or more 



rapidly using the Gribov-Migdal formula 29]. However, the contribution of the real part of 



/ into cross section is numerically small and thus can be neglected. 



B. Account of the running coupling constant. 



In order to obtain the full answer for the DGLAP approximation we also need to account 
for the radiative corrections to the propagators and vertices. An effective method to account 
for these corrections is to use the dispersion relations over the mass of the produced parton in 
the same way as in ref . 13|, |l6|, |30| . Our interest is in the kinematical domain: — t = A 2 > Ql 



otherwise dependence on t can be legitimately neglected. The integrand in eq. 13.51 has two 
potentially important kinematical regions: k\ < A 2 and k 2 3> A 2 . In the first kinematic 
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region the integrand is strongly suppressed and it is of order 1 / A 2 . The reason is that the IR 
singularities in the integrand are cancelled out when k t — > 0, k t — > A. As a result, the entire 
logarithmically enhanced contribution in the integrand comes from the second kinematic 
region. However, in this region 



2(kt,kt-A) 1 1 A 



(k t -A t )*(kt) {k t -A t f kf kf(A t -k 



Y2 




+ 7T2- T^T^T ~ Z-2 + O T4 ■ (3-9) 



The factor (kt, k t — A) follows from gauge invariance and accounts for the gluon polarizations 



271 ]. This means that within the logarithmic accuracy the integral over transverse momenta 



is the same as in the forward case, the only difference is that the integration starts at 



k 2 = —t. Hence we can directly use the results of ref. 16] . The renormalization of the ladder 



is given by diagrams of Fig. 3. Consider first the contribution of the exchanged gluon and 
two effective vertices. Taking the sum of discontinuities over two vertices and an exchanged 
gluon we obtain for the integral over transverse momenta 

(1/tt) f{d*k t /kl)dkllrn V{ ^ ^ * 7 . (3.10) 

Here as usual do, dp account for the multiplicative renormalization of the gluon propagator 

l/k 2 -> d G {k 2 ,a)/k 2 , (3.11) 

and renormalization of the propagator of the fermion 

l/k 2 -> d F {k 2 ,a)/k 2 . (3.12) 

(Above we omitted the additional terms proportional to 8{k 2 ) since they are not important 
for large t.) Therefore we can use the same approach as in ref. [l6|. Namely, we substitute 
k 2 — > —k 2 (l — z) and use kinematic identity: 

k 2 = -k 2 /(l-z) + zk 2 2 /(l-z). (3.13) 

Then the integral over k\ in the dispersion relation 13.101 is determined by the pole at k\ = 
—k 2 (l — z)/z. Thus in the Feynman integral for a renormalized cell we must put k 2 = 
-k 2 (l-z)/z. 



Since all renormalized vertices are explicitly written in ref. [16(, the only new element to 
take into account is that the left and the right vertices and the propagators dp depend now 
on different arguments. However since k 2 ^> A 2 integration does not introduce any changes 
relative to the case —t = 0. The factors a s (Q 2 ), dc, d F all combine together to the running 
coupling constant a s (k 2 ). Then the integral over transverse momenta in the ladder is given 
by 

d k t a s {k t )/k t = X = r^g{- — . x ■ ^ 2WA9 ), (3.14) 



log(-t + Qg)/A 2 ' 
where 

a ^ = ^SMWY b = u - 2Nllx (3 ' 15) 

We included Qq in eq. 13.141 to match parton distributions in the limit t — > 0. The same 
result can be obtained by calculating the integral in transverse momentum plane in the 
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whole region of k t explicitly, combining all three terms in the integrand together using the 
Feynman parametrization: 



, = /" d e ti A- {1 - x W +t{1 - x) \ (3.i 6 ) 

(k- (I- x)A) 2 - tx(l - x) 



o 

Finally let us note that we cannot distinguish in the LLA whether the argument in the 
coupling constant is k\ or k 2 . In the paper we follow suggestion of ref. [30], that choice of 
k 2 minimizes higher order corrections. 

The above results show that effectively the nonforward DGLAP distribution differs from 
the forward one only by the substitution of the variable \ = \ ^( log^/A 2 ) ) ^' gi ven by 
eq. 13. 141 Then the nonforward PDF satisfy the evolution equation: 

dG%(x, Ql Q\ t) _ a s (Q 2 - 1) _ ^ dz 



d\og(Q 2 -t) 4vr 



W -® B c {z)G c A {-,QlQ\t). (3.17) 
^ Jo z z 



Here &c(z) are the standard DGLAP kernels, G's are the nonforward ladders, A, B, C denote 
the parton species. Eq. 13.171 is valid for all massless partons A, B, C. The solution of this 
equation can be obtained from the solution for the forward distributions by the substitution 

X -» x'- 

The expression obtained in the previous subsection for the single cell in the approximation 
when the running of the coupling constant is neglected coincides with the expression in ref. 
This is because such DL terms are the same within DGLAP and BFKL approaches. 



13|| . In difference from ref. [26] we include the running coupling constant within the DGLAP 



framework and derive an evolution equation. 



C. Nonforward parton distributions for the hard inelastic exclusive diffractive 
processes within the LO DGLAP approximation. 

In order to evaluate the amplitude / we need to solve eq. 13.171 for the case of scattering 
off a parton. All we need to do is to solve the corresponding DGLAP evolution for zero 
t and then substitute y The corresponding parton distributions for forward case 



were derived in refs.[l3|, |16|, |3jJ: gluon distributions in the gluon, Dq(x, Q 2 ), and in quark 
Dq(x, Q 2 ). Since D q G <C Dq at small x [l6[ we shall need only the Dq function. 

As we already stressed above the pQCD evolution effects at small x are reduced to the 



double logarithmic(DL) approximation. Indeed, it was shown in ref. [15| that the PDFs 
evaluated using the DL approximation for gluon and quark distributions gives very good 
description of data at HERA energies. Moreover the double logarithmic approximation gives 



a good description of the parton distributions even at relatively large < 10~ 2 16|. As a 
result for the theoretical description of the small x inelastic diffraction it should be sufficient 
to use the double logarithmic expressions for parton distributions. In order to stress that 
we work with a gluon distribution with a single parton boundary conditions we shall denote 
this PDF as D(x,Q 2 ,t) below. 

The simplest description of the parton distribution in the DLA is achieved after the Mellin 
transform: 

D(j,t,Q 2 ,Q 2 )= f 1 z^ l D{x,t,Q 2 ,Ql). (3.18) 
Jo 
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The solution of the relevant DGLAP equation has the form 13j, [16 



D(j, X ) = exp( 1 (j) X ), (3.19) 

where 

AN 

iti) = —\ -a,a= (U/3)N C + (2/3)N f . (3.20) 

For our case N c = 3, Nf — 3. The solution in the x, x space is obtained after the inverse 
Mellin transform: 

xD(x,Q 2 ,t) = -L / djx- j G(j, X )- (3.21) 
Ztxi Jc 

The integration contour is chosen to the right of the j = 1 singularity and parallel to the 
imaginary axis. For the final answer we need to carry the inverse Mellin transform of eq. 13.191 
and to substitute x ~ > x' as ft was explained in the previous section: 

xD(x,Q 2 ,Q 2 ,t) =8N cX 'h(u)/u, (3.22) 

where 



u = ^lQN c \og(x/xj) X '. (3.23) 

The function I\ is the modified Bessel function of the first kind. For t — this is just 
the formula for Dq in the DLA, first obtained in ref. 3l|. For very small x we have the 
asymptotic expansion: 

/1/4 N l/A 

xD(x, x!) = —f=" r exp(u), (3.24) 

V A; V2¥log(xj/x) 3 / 4 PV ; V ; 

where xj is given by eq. 12.21 

The formulae 13.221 and 13.241 give the solution for nonforward parton distributions that 
will be building blocks for the diffraction amplitude /. 



D. The nonforward parton distributions beyond the DGLAP approximation. 



In the previous subsections we derived the formulae for the amplitude / within the 
DGLAP approximation (eq J3.22p which predicts the zero rate of the increase with energy 
for the nonforward PDF in the kinematical range — t > Q 2 + My. Formally we derived this 
result within LO DGLAP approximation. But account of NLO, NNLO, ...approximations 
will not change this prediction because all DL terms disappear in this kinematical domain. 

The DGLAP approximation ignores the possible contribution of ln(xo/x) terms that are 
not enhanced by large log(Q 2 /<5o — ^) terms. These terms are connected with the gluon 
radiation in the multi Regge kinematics. The BFKL and resummation models 
take into account both such terms and the double logarithmic terms 



17, 18, 32 



13L 117] . However, 



at present the direct numerical comparison between the BFKL and DGLAP results is not 
possible since the BFKL based models do not include the phase space constraints on the 
multi Regge gluon radiation, that follow from the energy-momentum conservation, and in 
addition neglect the running of the coupling constant. Recently the attempt was made to 



include such effects using the so called resummation models approach [18J, |22|, |32| . However, 
these were models developed for t = only. The interrelation between the Pomeron behavior 
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and the double logarithmic approximation in the multi Regge kinematics needs an additional 
investigation. 

In order to take into account the terms that do not contain large logarithms originating 
from the integration over parton transverse momenta it is convenient to represent a parton 
distribution as the product of DL contribution and the contribution that is not included in 
DLA. We may parametrize the nonforward parton distribution in the form 11.11 as: 

f(x, Q 2 , t) = f Uh (x, Q\ t)F(x, t). (3.25) 

As we already mentioned, the Pomeron behavior due to emission of gluon in the multi 
Regge kinematics is delayed till very high energies as a consequence of the strong constraints 
due to the local energy momentum conservation. 

Consider first the HERA energies. The interval in rapidity necessary for the emission of 



the gluons adjacent in rapidity in the multi Regge kinematics is at least ~ 2 4- 2.5 |33H35 
Since ]h(xq/x) > ln(My-/<5o), the photon fragmentation region occupies within the DLA at 
least > log(My/<5o) ^243 units in rapidity. The proton dissociation in the triple reggeon 
limit occupies at least three units in rapidity (due to acceptance of the HERA detectors). 
Rapidity span for the kinematics of HERA is ~ 8 4 9 units. Thus in the HERA kinematics 
it is hardly possible to emit even one additional gluon in the multi Regge kinematics. This 
means that at HERA energies if — t < Q 2 + M v there is no enough phase space for multi 
Regge corrections, and the resulted energy dependence is given the DGLAP terms, i.e. 
F = 1. For the kinematic range — t > Q 2 + M v the same kinematic analysis shows that 
because of the diminishing of the photon fragmentation region due to disappearance of DL 
terms where maybe a room for the radiation of one gluon within the multi Regge kinematics. 
We obtain: 

F(x, t) = l + p(Y- AY)6(Y - AY). (3.26) 

Here AY ~ 4 4 5 is the minimal rapidity needed for the start of a gluon radiation in a a 
multi Regge kinematics. Numerically we expect that the coefficient of the logarithmic term 
<C Qijp(t) — 1 since the existence of large logarithmic corrections of this type at HERA 
energies would contradict to a good agreement between DL asymptotics and experiment 



in the entire kinematic range of HERA [15j. The intercept of Pomeron trajectory ajp(t) is 
two - three times smaller in the resummation models, than in the LO BFKL approximation 
because of the accurate account of the DL terms 18, 32|. Thus with a good accuracy we 
can put F = 1 at HERA also for -t > Q 2 + M$. 

Consider now the energies significantly larger than the ones reached at HERA. In the 
ultraperipheral processes at the LHC one may reach x ~ 10 -6 , i.e. up to 10 units in rapidity 
available for the HID process, after the exclusion of the proton fragmentation region, that 
corresponds to 3 units in rapidity at least (see above). The good agreement between the 
results of DGLAP and resummation models up to x ~ 10~ 4 4 10~ 5 for Q 2 ~ 30 GeV 2 
(a s ~ 0.2) indicate that DL terms occupy the region of at least (2 4-3)log(Q 2 /<5o) in rapidity, 
i.e. 4 4 7 units for Q 2 ~ My ~ 10 GeV 2 for charmonium production. The DL terms define 
region in rapidity occupied by the photon fragmentation region and result in a reduction of 
the rapidity interval corresponding to the multi Regge kinematics. This means that only 
34-5 units in rapidity at most may be available for the emission of gluons in the multi Regge 
kinematics — less than the minimal region AY ~ 5. Thus in the ultraperipheral processes 
which could be studied at the LHC for small —t < My, our equations should be applicable, 
at least qualitatively. 
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For the kinematic range — t > Q 2 + M v the double logarithms are absent. In this case 
the rapidity range available for the gluon emission could reach 8-^9 units for multi Regge 
kinematics (after subtracting the photon fragmentation region (impact factor occupies ~ 2 
units in rapidity). Hence in this case 2-^3 gluons could be emitted in the multi Regge 
kinematics. 

Thus the DL behavior will continue to reveal itself in the entire interval of energies 
available at ultraperipheral collisions at the LHC Only at the maximum energies available 
at the LHC and —t > M v there is a sufficient phase space for the multi Regge gluon emission. 
In the kinematics of HERA the energy dependence is given by DGLAP approximation. For 
LHC energies, at the maximum energies to be probed in ultraperipheral processes, and at 
—t ~ My it will be interesting to look for the Pomeron behavior which will compete with 
the onset of the black disk regime. 

IV. THE ENERGY DEPENDENCE OF CROSS SECTIONS OF HARD INELAS- 
TIC PROCESSES WITH PROTON DISSOCIATION. 

Let us apply our results to the exclusive inelastic production of charmonium. The QCD 
factorization theorem 0, [1] allows us to evaluate the diffractive cross section in terms of the 
convolution of nonforward PDF discussed in the previous section with the parton distribution 
describing proton dissociation. The corresponding Feynman graph are shown at Fig. 1. We 
shall argue that the obtained above energy dependence of nonforward amplitude can be used 
as the interpolation formula correct also for — t ~ Q 2 + My. 

Let us consider first the kinematic range Ql < -t < Q 2 + My. Since -t < Q 2 + M v 
we may use the dipole approximation (tJ. We expand the impact factor into Taylor series 
over exchanged gluons transverse momenta , and then follow the steps used in the proof of 
QCD factorisation theorem. Thus the amplitude / of the process is proportional to 

f(-t,M v ,Q) = K J du J d 2 r(<f) V (u,r,A)A t (f) 1 (r,u,A)xD(x,4r 2 ,t). (4.1) 

The proportionality constant K is energy independent. It is matrix in the space of photon and 
vector meson polarizations: L — > L, T — > T,L — >■ T, T — >■ L, where L, T are longitudinally 
and transversely polarized quarkonium and incident virtual photon. 

D is a nonforward PDF, discussed in the previous sections. The argument of the nonfor- 
ward parton distribution x is given by eq. 12.101 

x ~ ( Xl + X2 )/2 ~ 3(M 2 + Q 2 )/(2xjs). (4.2) 

Let us stress that the only assumption used in the derivation of eq. 14.11 is the possibility to 
approximate impact factor by a first term in the Taylor expansion in exchanged gluon trans- 
verse momentum q^. Numerically such approximation is rather effective. The characteristic 
momenta of quarks in the impact factor r 2 ~ (Q 2 + My) /A (it is possible to neglect here 
the t-dependent term since -t C Q 2 + My). In addition u ~ 1/2 [36|. This means that the 
expression 14.21 can be actually approximated as 

/ = xD((M 2 + Q 2 )/(2x jS ), (My + Q 2 ), t)$(t, Q 2 , M 2 ), (4.3) 

where ^ 

$(t,Q 2 ,M 2 ) ~ J d 2 r J du((j) V (u,r,A)A t (j) 1 (r,u,A). (4.4) 
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$ depends on the quarkonium wave functions and influences the t-dependence of the cross 
section but it is independent of the incident energy. The energy dependence of the diffraction 
amplitude is entirely given by the convolution of nondiagonal forward distribution D((M V + 
Q 2 )/ (2xjs), M 2 + Q 2 , t) described in detail in the previous section with a parton distribution 
measured in DIS. Thus in the kinematical range Ql < -t < Q 2 +M 2 the energy dependence 
of the cross section at fixed Xj = —tjM\ is determined solely by the function D (eq. I3.22j) 
and does not depend on details of the photon and quarkonium wave functions. 

If momentum transfer — t ~ Q 2 + M v is large, the factorization theorem in the form 
described above should be modified. This is because the diffractive process is no more 
dominated by the DGLAP kinematics and the dipole approximation is not valid. The am- 
plitude is still given by the convolution of the impact factor and the gluon distribution. 
However in this kinematic region there is only one transverse scale — t, and the large loga- 
rithms originated due to QCD evolution are absent. Hence in the DGLAP approximation 
all amplitudes are reduced to energy independent impact factor. One would obtain the same 
result if we would simply extrapolate the energy dependence obtained above to this region. 
Thus we can use the energy dependence given by eqs. I3.22|I3.25I in the entire kinematic 
range Qq < — t < Q 2 + My. Beyond the DGLAP approximation the gluon distribution may 
contain the energy dependent terms which are not enhanced by large logarithms related to 
the Q 2 evolution. However it was argued in the previous section that such terms can give 
only small correction due to the energy-momentum constraints. 

We conclude that the energy dependence of the total diffractive cross section in the 
DGLAP approximation can be explicitly calculable in the kinematic range — t < Q 2 + M v 
as 

(h - $(t, Q\ My) 2 (4A ^ ))2 ((81/16)G(sj, t) + J2(qi(xj, t) + ft (xj, t)). (4.5) 



dtdxj " nil 



Here 



u = ^JlQN c \og(x/xj)x', X = t log 



1 , log(Q 2 /A 2 



b ^\og(-t + Q 2 )/A 2 " 
xj = -t/(M 2 x -m 2 p -t),x~ 3(Q 2 + My)/(2xjs), 6=11- 2/3N f , N c = 3, 

(4.6) 

and $(t, Q 2 , My) is the energy independent function, which depends on the details of the 
wave functions of the produced quarkonium and a photon. 

For the kinematic region -t ~ Q 2 + M v the cross section becomes energy indepen- 
dent in the kinematics covered by HERA, up to logarithmic corrections of the order 
(1 + 2f3(t) \og(x/x,j)) for \og(xj/x) > 4 5, and 0(t) < a F (t) - 1. 

The important case of charmonium photoproduction can be obtained from the previous 
formulae just by putting Q 2 = 0. 

To illustrate the pattern of the variation with t of the dependence of the cross section on 

i • i. • i • tt a i i ■ i • i • d\og(d 2 a I dtdx j) 

the rapidity gap interval we present m r lg. 4 the logarithmic derivative — : — - — for 

d\og(x/xj) 

the J/ijj photoproduction at —t = 0,4,9 GeV 2 . It corresponds to the effective "local" value 
of 2ajp(t) — 2. One can see that this quantity rapidly decreases with increase of —t. Note 
that the discussed approximation somewhat overestimates effective Q 2 for photoproduction 
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leading to a somewhat stronger energy dependence than a more realistic analyses of the 
t — photoproduction of J/ip j36j . 

Hence we conclude that for a fixed xj and — t > 4 GeV 2 the discussed process should 
depend very weakly on the rapidity gap interval at the HERA energies. This is consistent 
with the recent HERA findings [3| and phenomenological analysis of the earlier data 13 • 



It is worth mentioning that at sufficiently high energies where one may expect an increase 
of the amplitude due to the BFKL type dynamics the absorptive effects may strongly modify 
the energy dependence of the cross section. Really the process we consider is an example 
of the inelastic diffractive process. Such processes disappear in the black disk regime. In 
particular if xj is kept fixed and sufficiently large (say > 0.05) the contribution of the 
scattering at large impact parameters would be suppressed, while for the small impact 
parameters the probability of interaction with a gap in rapidity goes to zero. Since such 
effects strongly depend on the size of the QQ dipole they would also result in a weaker 
dependence of the cross section on t. These effects will be considered elsewhere. 

V. DVCS SCATTERING. 

The knowledge of the nonforward parton distribution evaluated above allows to calculate 
dependence on energy of the process of diffractive photon production: 7*(Q 2 ) + p — > 7 + p 
-DVCS. Although this process is hardly possible to observe at large t we present our results 
for completeness. 

Amplitude of this process is 

oc xG(x, Q 2 , t) = F(t, x, Q 2 ) xG(x, x'), (5.1) 

where xG(x, y) is the gluon distribution at HERA, whose energy dependence is well approx- 
imated by [15J 



xG(x, X ) ~ exp( v /l6A c log(x /x)x), x ~ 0.1, Q 2 ~ lGeV 2 , (5.2) 
and x' is given by eq. 13.141 

F 2g (t, x) = 1/(1 + t/m 2 (x, Q 2 )) 2 , (5.3) 
is dipole approximation for the two gluon form factor of the nucleon 37||. (We suppressed 



here the weak dependence of F 2g on Q 2 ). F 2g (t,x) may depend on energy in a restricted 
range of energies because of pion cloud around nucleon carrying small fraction of proton 
momentum. So its contribution at i > 0.1 should be negligible, But hard interaction 
rapidly increases with interaction energy and at sufficiently small x pion tail should reveal 
itself M . 



The expression 15 . 1 1 had been derived in the impulse approximation, which is questionable 
at large —t. However this effect is important for the calculation of absolute value of amplitude 
but it will not change the energy dependence. 

VI. CONCLUSIONS 

We found different QCD dynamics in the cross section of HID at HERA energies for fixed 
Xj = —t/Mx in the two kinematical regions. In the first region Qq < — t < Q 2 + My we 
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calculated the energy dependence of cross section within the DL approximation - eq. 14.51 
In the second region — t > Q 2 + M v we showed that the cross section of the HID processes 
is energy independent. This result is exact within the DGLAP approximation -valid within 
any order of LO, NLO....NNLO approximations. Our calculations explains observed in the 
recent HERA data on the HID processes [3| significant decrease of the intercept of pQCD 
"Pomeron" with increase of —t as compared to the intercept at t=0 . The corrections to 
DGLAP approximation due to gluon radiation in the multi-Regge kinematics are small at the 
kinematics of HERA, as it follows both from the analysis of energy- momentum constraints 



on the multi Regge gluon radiation, and the analysis of ref. |15J that DLA gives a good 
description of the HERA data. 

We argued that in ultraperipheral collisions at LHC the DL will give at least qualitatively 
good description of energy dependence of HID processes in most of the kinematic space 
available. However there may be a kinematic region x ~ 10~ 6 , -t > M v , that lies in the 
borderline of available energies where it will be very interesting to look for the onset of the 
pomeron behavior. 

In this paper we focused only on the energy dependence of the diffractive cross sections. 
The detailed numerical studies of the cross sections, including their absolute values and 



t-dependence will be presented elsewhere [39 



B.Blok is indebted to L.Lipatov for the illuminating discussion of the properties of BFKL 
and resummation models. 
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FIG. 2: The nondiagonal ladder 
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FIG. 4: The logarithmic derivative of the differential cross section as a function of energy for 
-t = 0,4,9 GeV 2 . 
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